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ABSTRACT 



We present new X-ray observations of the high mass X-ray binary (HMXRB) pulsar 
OAO 1657-415, obtained during one orbital period (10.44 days) with the Rossi X-Ray 
Timing Explorer (RXTE). Using the binary orbital parameters, obtained from Burst 
and Transient Source Experiment (BATSE) observations, we resolve the fluctuations 
in the pulse frequency at time scales on the order of one day for the first time. Recent 
BATSE results showed that OAO 1657-415 has spin-up/down trends in its pulse 
frequency time series, without any correlation with the X-luminosity at energies >20 
keV (Baykal 1997). In the present RXTE observations the source is found to be in an 
extended phase of spin-down. We also find a gradual increase in the X-ray luminosity 
which is correlated with a marginal spin-up episode. The marginal correlation between 
the gradual spin-up (or decrease in spin-down rate) and increase in X-ray luminosity 
suggests that the OAO 1657-415 is observed during a stable accretion episode where 
the prograde accretion disk is formed. 

Subject headings: accretion. X-ray binaries, OAO 1657-415 
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1. Introduction 

The high mass X-ray binary source HMXRB source OAO 1657-415 (OAO 1653-40) was first 
detected by the Copernicus sateUite (Polidan et al. 1978) in the 4-9 keV range. The HEAO-1 
observations also showed 38.22 sec pulsations in the 1-40 keV and 40-80 keV bands (White & 
Pravdo 1979, Byrne et al. 1981). Observations with Ginga and GRAN AT (Kamata et al. 1990, 
Gilfanov et al. 1991, Mereghetti et al. 1991, Sunyaev et al. 1991) have found pulse period 
changes. BATSE observations of this source with the Compton Gamma Ray Observatory (CGRO) 
showed that OAO 1657-415 is in an 10.44^^ binary orbit with an X-ray eclipse by the a stellar 
companion (Chakrabarty et al. 1993). The observed orbital parameters imply that the companion 
is a supergiant of spectral class B0-B6. The correlations between X-ray flux and pulse frequency 
derivatives (j>) fluctuations were investigated by using the previously published pulse frequencies 
and BATSE measurements (Baykal 1997). These correlations suggested that the formation of 
episodic accretion disks in the case of a stellar wind is the possible accretion mechanism. 

In this paper, we present the short term pulse frequency fluctuations and X-ray fluxes of OAO 
1657-415 in the light of recent RXTE observations. We have employed background subtraction by 
using the background models for the RXTE/PCA instrument and galactic ridge emission in the 
2-50 keV range. Our X-ray flux and pulse frequency measurements find an increase in the X-ray 
flux which is correlated with the decrease in the spin-down rate (or marginal spin- up trend) . 

2. Observation and Data Analysis 

OAO 1657-415 was observed between 1997 August 20-27 within the guest observer program 
of RXTE with proposal observation ID 20113. RXTE pointings of the source arc separated from 
each other by approximately six hours with a total observation span of 75 ksec. The results 
presented here are based on data collected with the Proportional Counter Array (PCA, Jahoda 
et al., 1996) and the High Energy X-ray Timing Experiment (HEXTE, Rothschild, et al., 1998). 
The PCA instrument consists of an array of 5 proportional counters operating in the 2-60 keV 
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energy range, with a total effective area of approximately 7000 cm^ and a field of view ~ 1° 
FWHM. The HEXTE instrument consists of two independent clusters of detectors, each cluster 
containing four NaI(Tl)/CsI(Na) phoswich scintillation counters sharing a common ~ 1° FWHM 
field of view. The field of view of each cluster is switched on and off source to provide background 
measurements. The net open area of the seven detectors is 1400 cm^ and each detector covers the 
energy range 15-250 keV. 

2.1. X-ray Light-Curves and Spectra 

Background light-curves and pulse height amplitudes are generated by using the background 
estimator models based upon the rate of very large events (VLE), spacecraft activation and 
cosmic X-ray emission with the standard PCA analysis tools. The background light-curves are 
subtracted from the source light-curve obtained from the Good Xenon event data (Fig. 1). Since 
the source is close to the galactic center, galactic ridge data were extracted from archival RXTE 
observations. Observations pointed at directions a few degrees away from the source are collected. 
The background spectra for the galactic ridge data were generated using PCA analysis tools. After 
the instrument background are removed from the galactic ridge data, the residual spectra are used 
as a background for the X-ray spectra of OAO 1657-415. The source spectrum was calculated 
using the same PCA background estimator models. (It should be noted that ~ 3 % systematic 
errors are used in spectral fitting of PCA data). For the HEXTE data the background subtraction 
is straightforward since the HEXTE detectors are rocking in 16 sec intervals. Standard ftools 
software for RXTE is used for the data reduction and for the dead time correction. We have used 
power-laws with high energy cut-off models together with a 6.65 kcV gaussian emission line which 
was found in previous GINGA observations (Kamata et al., 1990). We have found no evidence for 
any deviation from a power-law model with high energy cut off which might be attributable to a 
cyclotron feature. Table 1 presents the spectral parameters of X-ray spectra. Figure 2 presents the 
joint X-ray spectra of RXTE/PCA and HEXTE detectors. The power-law index, and the cut-off 
and e-folding energies are significantly different from those measured by Ginga, the powerlaw 
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steeper, but extending to higher energy. The broader energy band of the PCA and HEXTE 
combination better constrain the higher energy spectra, but the spectrum may vary intrinsically. 

2.2. Pulse Timing of OAO 1657-415 

The background subtracted light-curves are corrected with respect to the barycenter of the 
solar system. Using the binary orbital parameters of OAO 1657-415 from BATSE observations 
(Bildstcn at al., 1997), the light curves are also corrected for binary motion of OAO 1657-415 
(see Table 2). A long power spectrum was used to estimate the average pulse frequency. This 
pulse frequency is consistent with BATSE pulse frequency records at the same time (obtained 
through HEASARC http:cossc.gsfc.nasa.gov). In order to resolve the pulse arrival times and pulse 
frequencies at shorter time scales 29 pulse arrival times were generated (one pulse arrival time for 
each RXTE orbit). Pulse arrival times are found by folding the light-curve data into one average 
pulse for each RXTE orbit, folding all light-curves into one master pulse, and cross-correlating the 
master pulse with each of the 29 average pulses. In the pulse timing anaysis, we have used the 
method of harmonic representation of pulse profiles, as proposed by Deeter & Boynton (1985). 
In this method, pulse profiles are expressed in terms of harmonic series and cross-correlated with 
the master pulse profile. The maximum value of the cross-correlation is analytically well defined 
and does not depend on the phase binning of the pulses. The master pulse with 40 phase bins 
was represented by their harmonics (Deeter &; Boynton 1985) and cross-correlated with harmonic 
representations of pulse profiles from segments of the data. 

The pulse profiles for OAO 1657-415 axe variable. This affects the pulse timing. In order 
to estimate the errors in the arrival times, the light-curve of each RXTE orbit is divided into 
approximately 10-15 equal subsets and new arrival times are estimated. The average variance in 
the arrival times are computed and treated as errors of arrival times. The pulse arrival times are 
represented in Fig. 3. The residual pulse arrival times may arise from the change of the pulse 
frequency during the observation (or intrinsic pulse frequency derivative) and from the errors of 
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obital parameters (Deeter et al., 1981), 

,1., .o ^.asini 27r6T^/2 asini 27r asini , ^ , , 

dcp = (po+dl'{t-to) + -l'{t-to) +iyO{ )smln-V— COSln+V-^ dPorbit{tn-TT^/2)COSlr^ 

(1) 

where S4> is the pulse phase offset deduced from the pulse timing analysis, to is the mid-time of 
the observation, (p^ is the phase offset at to, Siy is the deviation from the mean pulse frequency 
(or additive correction to the pulse frequency), and z> is the pulse frequency derivative of the 
source, is the light travel time for projected semimajor axis, is the epoch when the 
mean orbital longitude is equal to 90 degrees, Pgrbit is the orbital period, S denotes the errors 
of these parameters, 1^ = 27r(f„ — / PorUt + is the mean orbital longitude at t„. The 
above expression is fitted to the pulse arrival times data. Table 2 presents the timing solution 
of OAO 1657-415 from our RXTE observations. The pulse frequency derivative obtained from 
the quadratic trend of the pulse timing analysis is vrxte = —(3.27 it 0.09) x 10""^^ Hz s~^. The 
average pulse frequency derivative we deduce from the pulse frequency history of BATSE archival 
data is — (3.1 ± 0.2) x 10~^^ Hz s~^, consistent at the la level. 

The pulse frequency records of OAO 1657-415 have shown that the source has stochastic 
spin- up/down trends (Baykal 1997, Bildsten et al. 1997) at timescales longer than weeks. In 
intervals between stochastic changes, OAO 1657-415 shows secular spin-up or down trends with 
lower values of noise strength (Baykal 1997). Our observations detect the source spinning down, 
on average. There are local deviations from quadratic trends which can be interpreted as very 
short term fluctuations (see Fig. 3). The torque noise analysis of OAO 1657-415 from the BATSE 
observations (Baykal 1997) showed that pulse frequency fluctuations at shorter than ~ 8 days are 
almost not detectable due to the measuremental noise, however RXTE observations are yielding 
significant fluctuations around the average quadratic trend (or spin-down rate), as shown in Fig. 
3. In RXTE observations we are able to construct pulses for time intervals as short as a few 
hundred seconds. This yields better timing at shorter time scales. Therefore the fluctuations in 
arrival times less than ~ 8 days are resolved. 
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2.3. Torque and X-ray luminosity changes of OAO 1657-415 

The X-ray flux and pulse frequency derivative correlations of OAO 1657-415 were investigated 
by using BATSE archival data base (Baykal 1997). BATSE pulse flux at 20-60 keV and pulse 
frequency series have shown no correlation between X-ray flux and pulse frequency derivatives. 
A strong correlation between specific angular momentum (1) and pulse frequency derivatives was 
found instead. These correlations implied that the specific angular momentum is directional, 
sometimes positive and sometimes negative (it 1), and that sometimes the flow is radial. These 
results suggested the formation of temporary accretion disks in the case of stellar wind accretion 
and the short term disk reversals are quite possible. In the analysis of BATSE data, the shortest 
time scales for resolving the significant pulse frequency fluctuations were of the order of 8 days. 

In the present work, in order to resolve pulse frequencies at shorter time scales, we used high 
counting statistics of RXTE/PCA detectors and obtained pulse arrival times at approximately six 
hours intervals. For each 4 or 5 pulse arrival times, we fitted a straight line segment and Tising the 
slope of this line we estimated a correction to the average pulse frequency {dcp = 4>o + ^^{t — tg))- 
In this way, we obtained pulse frequency records for each day roughly. In Fig. 4 we present 
the pulse frequency records estimated in this work. We accumulated the spectral data at 2-50 
keV, corresponding to observations of pulse arrival times. In each set of observations we fitted a 
power-law with high energy cut-off with a gaussian emission line, then we estimate unabsorbed flux 
(or intrinsic pulsar flux). Fig. 5 presents the pulsar flux history during the RXTE observations. 
During the observation, bolometric X-ray flux increased roughly by 70%. The spin-down rate is 
decreased and a marginal spin-up trend is seen. This is the flrst evidence in OAO 1657-415 to 
show positive marginal correlation between the X-ray flux and pulse frequency changes. 

3. Conclusion 

OAO 1657-415, has shown strong spin-up/down torques in its time history which can not 
be explained by wind accretion (Baykal 1997). The formation of temporary accretion disks was 
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therefore considered. If accretion onto the neutron star is from a Keplerian disk (Ghosh & Lamb 
1979), the torque on the neutron star is given by 

lu = n{ws)M Ik, (2) 

where Ik = {GMrg}^^^ is the specific angular momentum added by a Keplerian disk to the neutron 
star at the inner disk edge ro ~ O.Sr^ where = (2GM)~^/^/x^/^M~^/^ is the Alfven radius, /i is 
the neutron star magnetic moment, n{ws) ~ 1.4(1 — Ws/wc)/{l — Wg) is a dimensionless function 
that measures the variation of the accretion torque as estimated by the fastness parameter 
Ws = u/uK{ro) = 2'kP-^G-^''^M-^I'^ ^jp/'' M-^/'^ . Here Wc is the critical fastness parameter where 
the accretion torque is expected to vanish at Wc ~ 0.35 — 0.85 depending on the structure of the 
disk. According to this model the torque will cause a spin-up if the neutron star is rotating slowly 
{wg < Wc) in the same sense as the circulation in the disk, or spin-down, if it is rotating in 
the opposite sense. Even if the neutron star is rotating in the same sense as the disk flow, the 
torque will spin-down the neutron star if it is rotating too rapidly {wg » Wc)- In such a model 
one should see positive correlation between pulse frequency derivative (z>) and moderate mass 
accretion rate (M) if the disk is rotating in the same sense as the neutron star (Baykal 1997). 

Recent observations of accreting neutron stars have shown stochastic spin- up/down trends 
on time scales from days to a few years (Bildsten et al. 1997). Some of the sources switch 
from spin-up to spin-down states without showing great changes in their mass accretion rates 
(Bildsten et al. 1997). These unusual behaviors led Nelson et al. (1997) to the possibility of 
retrograde circulation of accretion disks. GX 4+1 shows correlation between the X-ray flux and 
the spin-down rate (Ghakrabarty et al. 1997) which may suggest a retrograde accretion disk. In 
our RXTE observations, GAG 1657-415 has shown marginal correlation with accretion rate and 
pulse frequency change. This positive correlation strongly suggesting that the disk formed in the 
spin-down episode is in prograde direction. Prom the BATSE observations, it was concluded that 
the pulse frequency derivatives and X-ray flux were not well correlated and the 8 days was the 
minimum for the correlation time scale (Baykal 1997). This RXTE observation implies that the 
time-scale of correlation is short, only a few days. To see the exact nature of correlations between 
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X-ray flux and pulse frequency derivatives, an even more extensive broad band X-ray observation 
should be carried out. 

A.B. thanks Ali Alpar and Mark Finger for critical reading of the manuscript and Jean .H. 
Swank, Tod .E. Strohmayer, Mike Stark, for stimulating discussions and USRA for supporting a 
visit to the GSFC. 
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Figure Caption 

Figure 1. The total(5 PCU) RXTE/PCA background subracted X-ray light-curve in 2-50 keV 
energy range. 

Figure 2. RXTE/PCA-HEXTE spectra of OAO 1657-415 (Note that HEXTE spectra is the 
summed spectra of Cluster 1 and Cluster 2). Below panel is the residuals of the fit in terms of 
values. 

Figure 3. Phase offsets in pulse arrival times. Solid line denotes the best fit of arrival times. 
Figure 4. Pulse frequency measurements of OAO 1657-415, from RXTE/PCA observations. 
Figure 5. 2-50 keV X-ray flux measurements of OAO 1657-415, from RXTE PCA observations. 
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Table 1: Spectral Fit Parameters for RXTE observations" 



Nh (102icm-2) 


12.94 ± 0.38 


Fe line center energy (keV) 


6.65 ± 0.03 


Fe line width (keV) 


0.40 ± 0.06 


Fe line intensity photons cm~^ sec ~^ 


(2.77 ± 0.31) X 10-3 


Photon Index 


1.07 it 0.03 


Cutoff Energy (keV) 


12.82 ± 0.29 


Folding Energy (keV) 


29.37 ± 0.98 


Power-law Norm. cm~^ at IkeV 


(5.96 ± 0.37) X 10-2 


Reduced 


1.28 (d.o.f. = 178) 



" Uncertainties in the spectral fit parameters denote single parameter 1 a errors. 



Table 2: Timing Solution of OAO 1657-415 for RXTE Observations" : 



Orbital Epoch (MJD) 


48515.99(5)'' 


Forb (days) 


10.44809(30)^ 


dkx sin i (It-sec) 


106.0(5)^ 


e 


0.104(5)'' 


w 


93(5)'' 


Epoch(MJD) 


50683.95400(2) 


Pulse Frequency (Hz) 


0.026775618(4) 


Pulse Preq. Derivative (Hz s"^) 


-3.27(9) X 10-12 



" Confidence intervals are quoted at the 1 a level. 

'' Orbital parameters are taken from Bildsten et al., (1997). Por6=orbital period, a^ sin i=projected 
semimajor axis, e=eccentricity, w=longitude of periastron. 
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